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ABSTRACT 


Calculated  reentry  aerodynamic  heating  effects  on  the  exter¬ 
nal  skin  of  an  ICBM  reentry  test  vehicle  having  a  relatively  low 
weight-to-dtag  ratio  arc  presented  in  this  report.  The  vehicle  is 
a  blunt  cone  with  the  aft  portion  designed  to  fail  just  after  maxi¬ 
mum  heating  and  the  forward  portion  designed  to  survive  to  impact. 

The  tape-wound,  reinforced  plastic  heat  shield  is  subjected 
to  maximum  reentry  heating  rates  between  125  x  10“*  and  650  x  104 
kcal/mzhr.  Resulting  external  and  internal  surface  temperature 
histories  are  given.  Effects  of  extreme  trajectory,  vehicle  char¬ 
acteristics,  and  atmosphere  variations  on  the  aluminum  substruc¬ 
ture  temperature  histories  are  discussed. 
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INTRODUCTION 

As  a  part  of  the  ARPA  experimental  reentry  radar  discrimination 
program,  the  U.S.  Army  Missile  Command  (USAMICOM)  has  been 
engaged  in  the  design  and  fabrication  of  several  typical,  lightweight 
ICBM  reentry  vehicles  having  various  ballistic  factors.  Vehicle  No.  1, 
which  is  discussed  in  this  report,  is  a  blunt  cone,  ablative  skin  vehicle 
designed  to  survive  intact  to  an  altitude  of  30.48  kilometers.  At  this 
altitude  the  forward  portion  of  the  vehicle  separates  xrom  the  destruct¬ 
ible  aft  frustum.  The  forward  portion  of  the  vehicle  is  designed  to 
survive  to  impact  so  that  recorded  data  can  be  obtained. 

The  purposes  oi  this  report  are: 

To  show  the  type  of  reentry  environment  to  which  the  vehicle 
is  exposed. 

To  show  the  calculated  ablative  heat  shield  requirements. 

To  predict  external  surface  temperatures  along  the  external 
surface  of  the  blunt  cone  vehicle. 

To  assist  in  determining  the  criteria  for  initiating  separation 
of  the  recovery  capsule. 

VEHICLE  DESCRIPTION  AND  ENVIRONMENT 

The  vehicle  is  a  blunt  cone  having  a  nose  radius  of  0.  15875  meter 
(6.  25  inches)  and  a  cone  half  angle  of  14  degrees  (Figure  1).  The 
initial  ballistic  coefficient,  W/CqA,  is  approximately  537  kg/m2 
(110  lb/ft2).  The  forward  portion  of  the  vehicle,  termed  the  recovery 
tip  in  Figure  1,  contains  a  tape  recorder  and  associated  recovery  gear 
and  is  designed  to  survive  reentry  to  impact.  The  aft  portion  of  the 
vehicle,  the  destructible  frustum,  is  designed  to  fail  at  a  reentry 
altitude  just  below  30.48  kilometers  (100,000  feet).  This  type  vehicle 
skin  design  is  advantageous  due  to  the  critical  weight  limitations  and 
other  vital  criteria  specified  in  the  vehicle  experiment  requirements. 

Shielding  of  the  reentry  vehicle  during  ascent  by  use  of  a  shroud 
is  considered  necessary  to  prevent  heat  from  soaking  through  the  heat 
shield  to  vital  components  before  the  reentry  phase  begins.  Nominal 
reentry  conditions  are  6,096  m/sec  (20,000  ft/sec)  at  an  angle  of 
23  degrees  below  the  local  horizon.  Shown  in  Figure  2  are  the  cal¬ 
culated  velocity  and  altitude  histories.  These  trajectory  parameters 
were  calculated  by  the  Advanced  Systems  Laboratory,  USAMICOM. 
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The  vehicle  spins  aboat  its  centerline  axis,  and  the  initial  angle-of- 
attack  is  less  than  5  degrees. 

THEORETICAL  PROCEDURES  AND  BASIC  ASSUMPTIONS 

Aerodynamic  heating  to  the  vehicle  was  calculated  by  real  gas 
procedures  based  on  the  widely  accepted  methods  of  Fay  and  Riddell, 
and  of  Rose,  Probstein,  and  Adams.  Aerodynamic  heating,  ablation, 
conduction,  pressure  distribution,  and  weight  calculation  procedures 
used  by  Stress  and  Thermodynamics  Analysis  Branch  are  described 
briefly  in  Appendixes  A.  B,  and  C.  To  obtain  a  simple  internal  config¬ 
uration  on  a  hemisphere  when  the  external  radius  and  required  thick¬ 
nesses  versus  angular  station  are  known,  a  method  was  developed 
whereby  the  weight  is  kept  to  a  minimum  while  the  specified  thick¬ 
nesses  are  met  or  slightly  exceeded  at  all  spherical  stations.  This 
procedure,  as  explained  in  Appendix  D,  is  for  an  internal  surface 
described  by  a  spherical  sector  with  the  radius  originating  on  the 
vehicle  centerline. 

Temperature  limits  were  set  on  the  aluminum  substructure,  as 
determined  by  the  sens. -ivity  of  internal  components  and  the  mechan¬ 
ical  properties  of  aluminum  alloys  versus  temperature.  For  the 
recovery  tip  the  design  temperature  is  approximately  378°  K,  and  for 
the  destructible  frustum  a  design  temperature  of  478°  K  was  selected. 
Based  on  these  substructure  limits,  parameter  studies  were  made  for 
each  vehicle  station  to  determine  the  thickness  of  ablation  material 
required  to  limit  the  aluminum  to  the  specified  temperature  limits. 

Ablation  calculations  are  based  on  the  assumption  that  the  fabri¬ 
cated  material  is  of  a  quality  at  least  equal  to  that  of  the  "lap-wound" 
glass -phenolic  produced  by  Westinghouse  Electric  Corporation  and 
tested  in  ABMA  aerodynamic  heating  simulation  facilities  in  1959 
and  I960.1'2 

The  initial  skin  temperature,  at  the  beginning  of  reentry,  was 
estimated  to  be  355°  K  to  account  for  heat  radiated  in  from  the  shroud 
during  ascent. 

The  net  effects  of  a  small  angle -of -attack  (5  degrees  or  less)  on 
the  heating  to  a  spinning  vehicle  are  assumed  to  be  negligible  so  that 
zero  angle -of-attack  procedures  can  be  used.  Some  investigators3 
state  that  this  approach  is  conservative. 
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HFAT  SHIELD  MATERIALS 


The  reinforced  plastic  heat  shield  selected  for  Vehicle  No.  1  is 
"lap-wound"  glass -phenolic.  This  material  was  selected  because  of 
its  low  cost,  availability,  suitability  to  tape  winding,  and  adequacy  as 
a  h  at  shield  in  a  relatively  mild  heating  environment. 

The  lap -winding  technique  is  extremely  flexible4  since  the  lami¬ 
nations  are  parallel  to  the  vehicle  centerline  and  a  flat  tape  is  used.  2 
The  orientation  of  the  tape  in  relation  to  the  gas  flow,  though  question¬ 
able  at  the  beginning  in  1958,  proved  to  be  of  no  importance  when 
tested  in  the  severe  environment  of  rocket  motor  exhaust  jets.  1,2 
The  ablation  results  for  lap-wound  materials  were  quite  comparable 
with  those  of  tape -wrapped  materials  having  laminations  making  a 
20 -degree  angle  with  the  cone  surface  (each  lamination  sloping  away 
from  the  centerline  toward  the  rear  of  the  vehicle). 

Glass -phenolic,  asbestos -phenolic.,  and  Thermolag  T-2305  were 
considered  for  protection  of  the  recovery  tip  buiinead  from  base 
aerodynamic  heating.  Of  these  materials  the  Thermolag  T-230, 
a  subliming  compound,  appears  to  be  the  most  practical,  weightwise. 
The  base  region  is  subjected  to  low  heating  rates  for  a  short  period 
of  time  after  separation  near  30.48  kilometers.  Since  the  shear  and 
heating  are  low  in  the  base  region,  a  reinforced  plastic  is  considered 
unnecessary  to  prevent  pitting  and  gouging. 

RESULTS  AND  DISCUSSION 

1.  Aerodynamic  Heating 


Calculated  real  gas  convective  heating  rates  to  a  hot  wall  are 
presented  in  Figure  3  for  the  recovery  tip  and  in  Figure  4  for  the 
destructible  frustum.  Maximum  reentry  aerodynamic  heating  occurs 
near  an  altitude  of  38  kilometers  (125,000  feet).  The  peak  heating 
rates  are  between  125  X  104  kcal/m2hr  on  the  aft  destructible  frustum 
and  650  X  104  kcal/m2hr  at  the  30 -degree  station  on  the  hemisphere. 
The  period  of  significant  reentry  heating  above  an  altitude  of 
30.48  kilometers  is  about  20  seconds.  Figure  3  shows  that  no  signifi¬ 
cant  aerodynamic  heating  occurs  after  45  seconds,  corresponding  to 
an  altitude  of  approximately  24  kilometers.  The  calculated  hea  fluxes 
on  the  recovery  tip  drop  quite  rapidly  after  38  seconds,  due  to  the 
sharp  decrease  in  velocity  resulting  from  a  significant  change  in 
ballistic  coefficient  when  the  recovery  tip  separates  from  the  aft 
frustum. 
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2.  Skin  Temperature  Histories 


Calculated  external  surface  and  internal  substructure  temper¬ 
atures  are  shown  in  Figures  5  through  14  for  several  vehicle  stations. 
For  the  recovery  tip  (Figures  5  through  10)  the  substructure  tempera¬ 
tures  peak  at  or  near  impact.  The  calculated  bulkhead  temperatures 
are  shown  in  Figure  11  for  Thermolag  T-230,  a  subliming  material  on 
the  outside  of  aluminum. 

In  Figures  12  through  14  are  calculated  external  and  internal  sur¬ 
face  temperature  histories  for  three  stations  on  the  aft  destructible 
frustum.  The  aluminum  substructure  temperature  rises  rapidly  near 
40  seconds,  due  to  the  recession  of  the  ablation  front  toward  the  sub¬ 
structure. 

Surface  temperature  distributions  at  several  altitudes  are  presented 
in  Figure  15.  Below  an  altitude  of  40.9  kilometers  (134,000  feet)  the 
entire  vehicle  surface  has  reached  the  ablation  temperature. 

3.  Ablation 

Ablation-depth  histories  calculated  at  selected  stations  along 
the  vehicle  are  shown  in  Figure  16.  As  expected,  the  forward  stations 
begin  ablating  at  altitudes  above  200,000  feet  and  the  aft  stations  begin 
ablating  at  altitudes  of  approximately  150,000  feet.  The  ablation  depths 
at  three  altitudes  are  plotted  versus  body  station  in  Figure  17.  From 
these  ablation  depths  the  weight  of  ablation  material  removed  versus 
time  or  altitude  was  determined  by  procedures  outlined  in  Appendixes 
B  and  C.  Figure  18  shows  the  calculated  ablation  weight-loss  versus 
flight  time.  At  100, 000  feet  the  total  weight  of  material  lost  is  almost 
25  pounds  or  15.  6  percent  of  the  initial  vehicle  weight.  Approximately 
5  pounds  have  been  removed  from  the  recovery  tip  and  20  pounds  from 
the  destructible  frustum  at  100,000  feet.  For  a  constant  CqA  this 
weight  loss  results  in  the  reduction  of  the  ballistic  coefficient  from 
110  lb/ft2  to  approximately  93  lb/ft2. 

4.  Heat  Shield  Requirements 

The  total  heat  shield  requirements  necessary  to  limit  the 
internal  substructure  to  378°  K  (220°  F)  on  the  recovery  tip  and  478°  K 
(400°  F)  at  100,000  feet  on  the  aft  destructible  frustum  are  shown  in 
Figure  19.  The  calculated  glass -phenolic  thicknesses  are  between 
14  and  24  millimeters  (0.  55  and  0.94  inch)  on  the  recovery  tip.  Thick¬ 
nesses  vary  from  3  millimeters  (0.  118  inch)  at  the  aft  end  to  4  milli¬ 
meters  (0.  153  inch)  at  the  fore  end  of  the  destructible  frustum. 
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Calculated  total  heat  shield  and  aluminum  substructure  weights 
are  tabulated  in  Table  I.  The  total  skin  weight  required  is  approxi¬ 
mately  65  pounds,  of  which  almost  45  pounds  are  on  the  destructible 
frustum.  Seventy  percent  of  the  total  skin  weight  is  attributed  to  the 
glass -phenolic  heat  shield. 

5.  Initiation  of  Separation  of  Recovery  Tip 


In  determining  a  criterion  for  initiation  of  separation,  several 
devices  such  as  timers,  temperature  sensors,  and  deceleration  switches 
were  considered.  The  deceleration  or  !,g"  switch  was  selected  as  being 
the  most  practical.  As  a  result  of  the  variation  of  "g"  levels,  aero¬ 
dynamic  loading,  and  subsequent  effects  of  aerodynamic  heating  on  the 
substructure  due  to  variations  in  trajectories  and  vehicle  parameters, 
additional  reentry  analyses  of  the  destructible  frustum  became 
necessary.  In  addition  to  the  nominal  trajectory,  minimum  and  maxi¬ 
mum  deceleration  flights  were  analyzed.  The  primary  differences  in 
the  trajectory  parameters  were  determined  by  the  Advanced  Systems 
Laboratory,  USAMICOM,  and  are  shown  in  Table  II. 

Calculated  aluminum  substructure  temperatures  versus  altitude 
for  the  destructible  frustum  are  shown  in  Figures  20,  21,  and  22.  As 
anticipated,  these  data  show  the  substructure  temperature  at  30.48  kilo¬ 
meters  to  be  highest  for  the  vehicle  flying  the  minimum  "g"  trajectory. 
The  substructure  temperature  is  lowest  at  any  given  altitude  for  the 
maximum  "g"  trajectory.  This  trend  is  reasonable  due  to  the  large 
variation  in  flight  time  from  100  kilometers  to  30.48  kilometers  and 
the  time  dependence  of  heat  flow  through  an  ablation  material  having 
a  low  thermal  diffusivity.  The  flight  time  from  100  kilometers  to 
30.48  kilometers  is  approximately  34.  5  and  25.  5  seconds,  respectively, 
for  the  minimum  and  maximum  "g"  trajectories. 

Immediately  below  an  altitude  of  30.48  kilometers  the  substructure 
temperatures  rise  sharply.  This  is  primarily  due  to  the  proximity  of 
the  receding  ablation  front  to  the  aluminum  substructure. 


CONCLUSIONS 

AMRAD  Vehicle  No.  1  is  exposed  to  a  relatively  mild  reentry 
environment.  The  major  portion  of  the  vehicle  is  exposed  to  maximum 
aerodynamic  heating  rates  of  300  X  104  kcal/m2hr  or  slightly  less. 
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The  calculated  glass -phenolic  heat  shield  required  for  Vehicle 
No.  1  weighs  about  46  pounds  with  62  percent  of  this  located  on  the 
destructible  frustum.  The  aluminum  substructure  weighs  approxi¬ 
mately  19  pounds  with  87.  5  percent  of  this  located  on  the  destructible 
frustum. 

At  any  given  altitude  above  30  kilometers,  the  calculated  sub¬ 
structure  temperature  for  the  aft  destructible  frustum  increases  as 
the  vehicle  trajectory  goes  from  the  maximum,  to  nominal,  to  minimum 
deceleration  cases.  The  increase  in  flight  time  to  a  given  altitude  is 
the  major  factor  in  this  trend. 


Table  I.  Skin  Weights 


Vehicle  section 

Aluminum 
thickness 
(in.  ) 

Aluminum 

weight 

<ib) 

Heat  shield 
weight 
(ib) 

Total  skin 
weight 
<ib) 

Recovery  tip 

0.054 

2.485 

17.423 

*19.908 

Destructible  frustum 

0.  054 

16.  62 

28.  28 

44.90 

^Includes  insulation  and  substructure  on  recovery  tip  bulkhead. 


Table  II.  Relationships  Between  Maximum  and  Minimum 
Deceleration  Trajectories  Compared  With  the  Nominal 


Parameter 

Maximum  "g" 

Minimum  "g" 

w/cda 

90  %  nominal 

110%  nominal 

Reentry  velocity 

104%  nominal 

95%  nominal 

Reentry  angle 

26  degrees 

20  degrees 

Wind 

Headwind 

Tailwind 

Recovery  tip* 


■  Destructible  frustum 


0.8666  m 


(K33  mJIdO3  ft) 


Figure  2.  Reentry  Trajectory  for  Vehicle  to  30. 48  Kilometers 
and  for  Recovery  Tip  to  Impact 


Calculated  convective 


Figure  5.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Stagnation  Point  of  the  Recovery  Tip 


V.  «  6096  m/»*c 
6  *  23° 


Figure  6.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  30 -Degree  Station  of  the  Recovery  Tip 


Figure  7.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  55 -Degree  Station  of  the  Recovery  Tip 


Figure  8.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  76 -Degree  Station  of  the  Recovery  Tip 


Ablotlnq  temperature 


V-  *  6096  m/«eo 

e\ * 23® 

Slot  Ion  L*  «  0.236 


0  L*  =  0236  m 


-R  s  013875  m 


300  100 


Gloee-phenollc  thlckneee  *  14.4  mm 
Aluminum  thlckneee  *  1.372  mm 


0  80  120  160 
Reentry  flight  time  from  122.0  km,t  (tec) 


-Altitude  (M t) 


Figure  9.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Lx  =  0.  236 -M  Station  of  the  Recovery  Tip 


Ablotlnq  temperature 


V#  *  6096  m/eec 
8,  -■  23° 

Stotlon  '  Lx  «  0,33 


3  L,  s  0.35  m 


■  R  *  013073  m 


Gl o« i -phenolic  thlckneee  =  14.29  mm 
Aluminum  thlckneee  *  1.372  mm 


-Altitude  (kft) 


Reentry  flight  time  from  122,0  km,  t  (tec) 


Figure  10.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Lx  =  0.  35 -M  Station  of  the  Recovery  Tip 


Figure  11.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  on  the  Bulkhead  of  the  Recovery  Tip 


Figure  12.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Lx  =  0.4-M  Station  of  the  Destructible  Frustum 
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Figure  13.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Lx  =  0.  75-M  Station  of  the  Destructible  Frustum 


Figure  14.  Calculated  Exposed  and  Unexposed  Surface  Temperature 
Histories  at  the  Lx  =  1.  4-M  Station  of  the  Destructible  Frustum 


Calculated  weiglit  of  ablation  material 


V„  «  6096  m/ftcc 
6,  =  83* 

Mot*rid:  "lop -wound"  gu»i-ptono!ic 


Figure  17.  Calculated  Ablation  Depths  at  Three  Altitudes 

Versus  Body  Station 


Figure  18.  Calculated  History  of  Weight  Removed  by  Ablatio 


Figure  19.  Calculated  Total  Insulation  Thickness  Requirements 


Figure  20.  Calculated  Substructure  Temperature  Histories 
for  Forward  Station  of  Aft  Frustum 


ucture  temperature 


Figure  22.  Calculated  Substructure  Temperature  Histories 
for  Aft  Station  of  Frustum 
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Appendix  A 

USAMICOM  AERODYNAMIC  HEATING  METHODS 

Widely  accepted  real  gas  methods  are  employed  in  the  USAMICOM 
theoretical  heating  analyses.  At  the  stagnation  point  the  Fay  and 
Riddell  method  with  a  modified  Lewis  number  (L)  contribution6  is  used. 

The  equation  for  qc  at  $  =  0°  is 

°.  i  0.4/  v  0.25T  .  _h  vl0.48 

<lc=CiDn  (Pw^w)  (P^g)  |  1+(L“1)^  ~——J  (hs-hw) 

where  Cx  is  a  constant. 

At  body  stations  other  than  the  stagnation  point,  the  laminar  heating 
rate  is  found  by  modifying  the  stagnation  point  rate  according  to  the 
pressure  distribution. 


For  turbulent  flow  the  theoretical  heating  rates  are  based  on 
r:  ethods  developed  by  AVCO7  with  the  modified  L  contribution  inserted. 


To  be  conservative  in  the  design  of  hardware,  both  laminar  and 
turbulent  heating  rates  are  calculated  at  any  given  station  except  the 
stagnation  point.  The  higher  heating  rate  is  then  selected  for  use  in 
determining  the  structure  response.  In  this  manner,  transition  is 
effected  and  no  other  margins  of  safety  are  incorporated. 

Ideal  gas  procedures  based  on  methods  of  Sibulkin,  Eckert,  and 
Van  Driest  have  also  been  used  extensively  in  the  past.  By  proper 
modifications  the  reference  temperature  ideal  gas  procedure  gives 
heating  rates  comparable  to  the  real  gas  procedure.  8 

The  pressure  distribution  on  a  hemisphere  is  determined  from 
the  modified  Newtonian  impact  theory  matched  with  a  Prandtl-Meyer 
expansion  to  the  tangent'point.  For  the  frustum,  a  blast  wave  method 
is  used  to  find  the  pressure  decay  between  the  tangent  point  and  the 
sharp  cone  value.  These  procedures  are  slightly  conservative 
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compared  4o  the  GASL  3-D  method  of  characteristics. 9  Results  from 
the  blunt  cone  pressure  routines  used  herein  have  been  favorably  com¬ 
pared  with  experimental  pressure  data.9 

The  viscosity,  compressibility  factor,  specific  heat  ratios,  local 
enthalpy,  local  temperature,  local  speed  of  sound,  and  local  density 
are  taken  from  data  published  by  C.  G.  Hansen10  and  J.  Hilsenrath. 11 

On  the  hemisphere  a  modified  Newtonian  flow  theory  is  used  to 
determine  the  velocity  gradient. 

Before  and  after  ablation,  a  forward  finite  difference  heat  conduc¬ 
tion  procedure  is  used.  During  ablation  a  special  conduction  procedure 
applicable  to  conduction  in  a  material  having  a  receding  surface  is 
utilized.  Both  of  these  procedures  compare  well  with  exact  solutions. 

Ablation  of  material  in  flight  is  based  on  experimental  results  of 
material  performance  tests  in  liquid  propellant  rocket  motors.  The 
average  effective  heat  of  ablation,  H^»  versus  heating  rates  (Figure  23) 
is  determined  by  use  of  ablation  depth  measurements  and  metallic 
calorimeters. 


Heat  flux  to  nonablating  surface  at  the  elating  temperature 

Figure  23.  Experimental  Thermal  Performance 
of  a  Reinforced  Plastic 
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Ablation  and  calibration  models  used  in  the  tests  had  tip  diameters 
from  19.05  to  635  millimeters.  The  stagnation  temperatures  were 
from  1920°  to  4700°  K  and  the  pressures  were  between  70, 000  and 
105,000  kg/m2.  Hot  wall  heat  fluxes  ranged  between  50  X  104  and 
2000  X  104  kcal/m2hr. 

Even  though  the  rocket  motor  exhaust  gas  composition  differs  con¬ 
siderably  from  air,  excellent  flight  and  theoretical  ablation  correlations 
have  been  obtained  using  this  semiempirical  method.  One  subscale  and 
three  full-scale  IRBM  vehicles  flying  near  Mach  15  were  recovered 
after  flight  and  measured  for  material  removal  during  the  JUPITER 
research  and  development  program.  12' 13' 14  These  ablation  depth  corre¬ 
lations  were  for  glass -melamine,  glass -phenolic,  refrasil-phenolic, 
and  asbestos -phenolic  materials.  The  flight  ablation  depths  of  all  four 
materials  agreed  very  closely  with  the  predicted  values,  especially  on 
the  hemispherical  portion  of  the  blunt  cones. 

Good  comparisons  have  also  been  obtained  between  theoretical  and 
flight  data  for  the  heat  shield  on  a  BIG  JOE  capsule  reentering  the 
earth's  atmosphere  approximately  6,096  m/sec.14'15 
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Appendix  B 

VOLUME  OF  HEMISPHERICAL  SECTION 


Often  in  heat  transfer  work  and  in  the  design  of  blunt  vehicles,  it 
is  necessary  to  calculate  the  weight  of  material  removed  by  ablation 
or  to  determine  the  entire  skin  weight  of  hemispherical  tips.  To 
obtain  the  weight  it  is  required  that  the  volume  be  known. 


The  volume  of  the  spherical  section  shown  in  Figure  24  is  readily 
computed  with  the  following  equation: 

n  =  4/Afl 

Vt  =  ZttAQ  l  [R2Q6n  -  R06?,  +  6^/3]  sin  6n  (1) 

n=  1 


Figure  24.  Volume  of  a  Spherical  Section 


This  equation  is  derived  as  follows: 

x  is  the  radius  vector  to  the  center  of  gravity  of  the  incremental 
area 
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y  =  x  sin  Q  is  the  moment  arm 
dA  =  xdfidx  is  the  incremental  area 
dV  =  2ffydA  =  2nxz  sin  0d@dx 


Jro  =  4>  /*rq 

I  x£  sin  0d0dx 

9  =  0°JRn  - 


=t  r  *  k~  (ro  -  sinede 

*  Je  =  o° 


v*  = 


n  =  4>/A0  r  , 

ZnAe  l  [r206„  -  R04.  +  63n/3j  sin  6n 
n  -  1 


(2) 

0) 


By  proper  selection  of  the  incremental  angle,  A0,  commensurate 
with  the  particular  variation  of  wall  thickness,  accurate  volumes  on  a 
spherical  tip  can  be  obtained  by  use  of  Equation  (3).  The  thickness 
6n,  measured  normal  to  the  external  surface,  is  the  average  thickness 
over  the  A0  at  increment  n. 
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Appendix  C 

VOLUME,  CALCULATION  FOR  FRUSTUMS 
OF  RIGHT  CIRCULAR  CONES 

Equations  determining  the  volume  for  frustums  of  right  circular 
cones  were  derived  using  Pappus'  second  proposition: 

"If  a  figure  of  area  A  revolves  about  an  axis  in  its  plane  but 
not  cutting  it,  then  for  a  complete  revolution,  the  volume  of  the  solid 
generated  is: 

V  =  ZnyA  (4) 

where  y  is  the  distance  from  the  axis  to  the  center  of  gravity  of  A.  " 

The  frustum  volume  equations  were  deduced  for  a  linear  and  a 
nonlinear  variation  of  skin  thickness,  6,  with  length,  Lx. 

1.  Linear  Variation  of  Thickness  With  Length 

The  volume  of  the  frustum  of  a  right  circular  cone  having  a 
linear  variation  of  thickness  with  length,  Lx,  may  be  calculated  by 
adding  the  volume  of  the  solid  generated  by  the  rectangular  area  "A" 
(Figure  25)  to  the  volume  of  the  solid  generated  by  the  triangular  area 
"B"  in  the  following  manner: 

For  area  "A"  the  volume  is 

Lx  <5, 

VA  =  27t62Lx  R  sin  Q  +  —  cos  0  -  -Z-  sin  Q  (5) 

and  for  area  "B"  the  volume  is 


Combining  Equations  (5)  and  (6)  gives  the  total  volume,  Vt 
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Figure  25.  Linear  Variation  of  Thickness  With  Length  (L) 


2.  Nonlinear  Variation  of  Thickness  With  Length 


The  volume  of  the  frustum  of  a  right  circular  cone  having  a 
nonlinear  variation  of  thickness  with  length  can  be  calculated  from  the 
following  equation  where  n  is  the  increment  being  considered  in  the 
x  direction: 


n  =  x^/Ax 
Vt  -  ZttAx  £ 
n  =  1 


R6n  + 


cos  @Xji  6n 
sin2  B 


Sn 


(8) 
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In  deriving  Equation  (8)  (Figure  26): 


The  incremental  area  is  dA  =  6dL 
The  moment  arm  is  yc 

The  incremental  volume  is  dV  =  27ryc6dL 


Figure  26.  Nonlinear  Variation  of  Thickness  With  Length  (L) 


Other  pertinent  equations  are: 


dL 


dx 

sin  0 


Yj  =  R  sin  0 
Yz  =  Yl  +  x  cot  Q 
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Yc  =  Y2  -  j  sin  0 

Y  =  R  sin  0  +  x  cot  6  -  y  sin  0 

Equation  (8)  can  be  used  to  accurately  determine  the  volume  of 
ablation  material  removed  due  to  hypersonic  flight  in  the  atmosphere 
and  the  volume  of  an  entire  frustum  of  nonlinear  thickness  with  length. 
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Appendix  D 


SIMPLE  INTERNAL  SURFACE  DESCRIPTION  FOR  A  HEMISPHERICAL 
TIP  HAVING  A  SKIN  THICKNESS  VARYING  WITH  ANGULAR  STATION 

In  the  design  of  reentry  vehicles,  it  is  frequently  desirable  or 
necessary  to  deviate  slightly  from  optimum  ablation  material  or  heat¬ 
sink  thickness  requirements  on  a  hemispherical  tip  to  simplify  the 
internal  surface  configuration.  A  simple  internal  configuration  is 
advantageous  due  to  reduced  costs  of  fabrication,  especially  for  ve¬ 
hicles  that  are  to  be  produced  only  in  limited  numbers. 

For  a  given  external  spherical  configuration  and  for  known  skin 
thickness  requirements  varying  with  angular  stations  (Figure  27),  a 
procedure  has  been  derived  for  determining  an  internal  surface  con¬ 
figuration  described  by  a  spherical  sector  with  the  radius  originating 
on  the  vehicle  centerline. 


Figure  27.  Procedure  for  Describing  a  Simple  Internal 
Surface  Configuration 


Known  parameters  required  for  solution  of  the  internal  configura¬ 
tion  are  (Figure  27): 
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R0  -  External  radius 

-  Skin  thickness  measured  normal  tc  the  external  surface 
at  any  known  0j 

-  Skin  thickness  measured  normal  to  the  external  surface 
at  any  known  02 

It  is  required  to  find: 

Center  location,  C,  on  the  centerline  for  the  internal  surface 
Internal  surface  radius, 

Calculated  skin  thicknesses  at  any  desired  angular  station 

With  the  known  parameters  listed  in  Figure  27  let 

Xj  =  (Ro  “  61)  cos  03 
X2  =  (Ro  “  62)  cos  0 2 
Yi  =  (Ro  “  61)  sin  0j 
Yz  =  (Ro  “  62)  sin  02 


Then  from  triangle  ACD 


Yi  +  (xl  +  Ax)2  =  R- 

(9) 

and  from  triangle  BCE 

Yz  +  (x2  +  Ax)2  =  Ri 

(10) 

Solving  Equations  (9)  and  (10)  simultaneously  results  in 
Y\  -  Yz  +  (*3  +  Ax)2  -  (x2  +  Ax)2  =  0 


or 

(xi  +  y2)  -  (xf  +  yf)  +  2Ax(xj  -  x2)  =  0 

and 

Ax  =  (x2  +  yi)  ~  (Xi  +  Yi) 

2  (xj  -  x2 ) 


(11) 


(12) 
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Equation  (12)  defines  the  location  of  the  internal  surface  center  lying 
on  the  vehicle  centerline. 

Substituting  the  Ax  from  Equation  (12)  into  Equation  (9)  yields  the 
solution  for  the  internal  radius 

Ri  =  yjYi  +  (xi  +  Ax)2  (13) 

To  find  the  skin  thickness  at  any  spherical  station,  0 ,  let 
x  =  (R0  -  6q)  cos  0 

and 


y  =  (R0  “  fig )  sin  0 

then  from  right  triangle  relations 

y2  +  (x  +  Ax)2  =  R?  (14) 

where  Ax  is  length  found  from  solution  of  Equation  (12).  Substituting 
x  and  y  into  Equation  (14)  gives 

(Ro  ~  60 )2  sin2  0  +  [(R0  “  60 )  cos  0  +  Ax]2  =  R? 

which  expands  to 

(Ro  "  60 )2  +  2Ax(Rq  -  60 )  cos  0  +  (Ax2  -  R?)  =  0  (15) 

Solving  for  RQ  -  6q  in  Equation  (15)  yields 

-2Ax  cos  0  ±  vAax2  cos2  0-4  ^x2  -  Ri) 

K0  -  60  =  - - -  (16) 


Since  RQ  -  6q  must  be  positive  for  the  hemispheres  considered,  the 
radical  in  Equation  (16)  must  be  positive;  thus  solving  for  5q,  the 
skin  thickness,  at  any  angular  station,  0,  gives 


60 


R0  +  Ax  cos 


sin2  0 


(17) 


To  quickly  obtain  the  best  solution  of  an  internal  hemispherical 
configuration  with  the  center  on  the  vehicle  centerline,  several  pairs 
of  required  skinthicknesses  (determined  by  stress  and  thermodynamic 
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analyses)  at  known  angular  stations  are  investigated  by  use  of  a  high¬ 
speed  digital  computer.  The  calculated  skin  thickness  obtained  from 
the  computer  program  for  each  internal  surface  configuration  (one 
configuration  for  each  pair  of  skin  thicknesses  input  to  the  program) 
is  compared  with  known  skin  thickness  requirements  previously  deter 
mined  from  stress  and  thermodynamic  analyses. 

Selection  of  the  best  internal  surface  configuration  is  then  made 
from  the  skin  thickness  comparisons  by  considering  primarily  the 
effects  of  skin  thickness  deviations  on  accomplishment  of  the  mission 
and  the  effects  of  skin  thickness  deviations  on  the  total  weight  of  the 
hemispherical  section. 

With  the  skin  thickness  incorporated  as  an  output  to  the  configura 
tion  procedure,  the  weight  of  the  material  on  a  hemisphere  is  readily 
obtained  by  use  of  the  volume  equations  presented  in  Appendix  B. 
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